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Time ordering may be defined by first defining the limit of no time ordering (NTO) in terms of a 
time average of an external interaction, V(t). Previously, time correlation was defined in terms of 
a similar limit called the independent time approximation (ITA). Experimental evidence for time 
correlation has not yet been distinguished from experimental evidence for time ordering. 
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I. INTRODUCTION 

A. Space and time 

Both space and time are parameters used to mathematically describe observable properties of physical objects or 
systems of objects. In some ways space and time are similar, both intuitively and mathematically since they are the 
most bàsic coordinates used to describe the dynamics of physical systems. It is oftcn convenient and conventional to 
use t and r to specify when and where an object is, even in quantum systems where the precision (or locality) of these 
coordinates is limitcd by the unccrtainty principle or obscured by entanglement. Nevertheless, space and time diffcr 
fundamentally (even neglecting the second law of thermodynamics, which dehnes a direction in time but no preferred 
direction or arrow in space). Because there is no counterpart of temporal causality in a space-like context, an object 
can repeatedly return to any spatial location, but it cannot return or jump ahead temporarily in time. Much has been 
written about spatial properties of multi-particle systems, but less about time, even though thrcc-dimensional space 
would seem much harder to deal with than one-dimcnsional time. Time, unlike space, seems somehow enigmàtic. 
Newton refers to time as "like a river flowing." Einstein calls time "that which a clock measures." And Feynman refers 

to time as the space between events, i.e. that which keeps EVERYthingFROMhappeningALLatoNCE. 

One purpose of this paper is to see what happens when we try to carry an idea formulated in a spatial context 
into a time context. That idea is correlation (associated with non-randomness, or quantum entanglement). Time 
correlation 0, Q has been directly related to time ordering, a causality condition implied by the time-dependent 
Schròdinger equation that constrains sequential interactions to occur in order of increasing time. 

B. N-body problem 

It is well known that the N-body problem, namely determining the quantum evolution of N correlated or interacting 
bodies, is exponentially difHcult. Kohn |3( specifically estimates that the amount of computer storage capacity needed 
to solve an N-body problem scales as e 3N , where 3 is simply a coeficient based on experience. Combining Kohn's 
estimate with Moore's law, which estimates, again based on experience, a doubling of computer capacity every 1.5 
years, one can easily show that 6.5 years are required for computer capacity to grow enough to accommodate one 
additional body. Thus, increasing capacity from, say helium (N = 3 ignoring nuclear and subnuclear structure) to 
carbon (N = 7) is estimated to require about 25 years, progressing to wàter (N = 11 ignoring the inner shell) requires 
about 50 years, and DNA thousands of times the age of the universe. The N-body problem is difficult. On the 
other hand, since much of physics, chemistry, biology, materials science, nanotechnology, and quantum computing 
involves múltiple partides, the problem is significant. In quantum computing, for example, much has been understood 
about individual qubits (regarded here as single "bodies"), but much less is understood about the coupled networks 
of qubits needed to build a quantum computer. Evidently, there is a need for sensible, well-defined approximations 
to the N-body problem. Perhaps the most widely used of these is the independent particle approximation (IPA), 
where the N-body problem is dramatically reduced in complexity to a problem of N independent bodies (or partides, 
quasi-particles, qubits, . . . ). 

C. Correlation 

It is sometimes argued that a theory is not well defined until a reliable method is developed to calculate corrections 
to the theory. In the case of the IPA, the corrections are called correlations and represent interconnections among 
the positions of the various partides. Correlation implies complexity: a system of independent (i.e. uncorrelated) 
partides is less complex than a system in which each particle's behavior depends on the behavior of every other 
particle. From one viewpoint, correlation is a key to understanding how to make complex systems from simple ones. 
From the opposite point of view, in a very complicated system correlation can be a pathway to see order through a 
landscape of chaos. 

Our approach is to focus on the first viewpoint - building complex systems from simple ones. If the correlations are 
sufhciently small, then the IPA is reasonably accurate. Correlation, one method to approach the N-body problem, 
is defined as the difference between a full solution and the IPA. Thus, the IPA limit of no correlation is used to 
define correlation itself. While this approach appears a bit awkward at first, it is conventional in the study of non- 
random processes and quantum entanglement as well. For example, in the limit of no entanglement the N-body wave 
function may be written as a simple product of independent particle solutions, and entanglement may be defined as 
the deviation from this limit. 
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In most applications, correlations arise from spatial interactions between partides. In N-electron àtoms, for example, 
electron correlation is generated by the l/r« Coulomb interactions between electrons. In this context, the IPA is also 
known as the independent electron approximation (IEA) and is defined by replacing these complicating correlation 
interactions by a sum of single-electron mean-field interactions. The resulting Hamiltonian is a sum of single-electron 
terms, even in the case of dynamic atòmic systems jj], and the resulting wave function is a product of single-electron 
wave functions. One common method for analyzing deviations from the IEA is many-body perturbation theory 
(MBPT), used here and in other fields as well to describe correlation at various orders in a given perturbation. 

One of the primary qüestions motivating the work described here is, "what happens if one tries to define an indepen- 
dent time approximation (ITA) , where one eliminates correlations in time rather than the conventional correlations in 
space?" In such an approximation, múltiple times may be used for múltiple partides 0- This is simpler than using a 
single time for all partides, just as the use of N independent particle positions is simpler that trying to solve a single, 
but complicated, N-body problem. Time correlation may then be defined as a deviation from the ITA in analogy to 
the dcfinition of spatial correlation by reference to the IPA. 

D. Correlation in time and time ordering 

In this paper we first review early efforts to define the ITA, in which a key step is to remove time ordering between 
fields acting on different partides. This is closely related to the approximation of no time ordering (NTO), where 
all time ordering is removed for fields acting either on different partides or on the same particle. These early efforts 
to define and study the ITA were based on a perturbation expansion of the external interaction V(t). Since time 
correlation or time ordering first appear at second order in V(t), the observable effects are small and the calculations 
difficult. Next we review more recent investigations of time ordering and the NTO approximation for qubits in 
non-perturbative external fields. Therefore, the organization of this paper is unusual. The harder problem of time 
correlation is addressed first, and the related but more tractable problem of time ordering is discussed second. We 
follow this unusual order because it reflects what has happened in recent years and also because the newer work on 
time ordering raises some intriguing qüestions and challenges. In particular we emphasize that experimental evidence 
for time correlation has not yet been distinguished from evidence for time ordering. That is, no difference between 
the ITA and the NTO approximations has yet been observed in any experiment, although recent experiments have 
been progressing in this direction, as we illustrate below. 

II. THE N-BODY PROBLEM 

In this section we consider solutions to the N-body problem described by a Hamiltonian 

H = Ho + Vs(t), (1) 

where there may be many terms in both the unperturbed Hq, which we assume to be solvable, and in Vs(t), which 
we regard as an external time-dependent interaction (written in the Schròdinger representation) . Where appropriate, 
Vs(t) may be treated as a classical external field. In the case of atòmic collisions, this corresponds to the widely used 
semiclassical approximation (SCA) where the projectile is regarded as a classical particle if it is a proton or electron 
and a classical wave if it is a photon. 

The N-body problem may be solved in various representations, depending on how H is separated into Hq and Vs(t), 
ranging from the Schròdinger representation where effectively Vs (i) — H and Hq = to the Heisenberg representation 
where Ho = H and Vs(t) — 0. In the SCA, singularities can occur Q in the first order amplitudes when the Schròdinger 
representation is used. Except where otherwise specified we shall use an appropriate intermediate representation that 
takes maximum advantage of known solutions for Ho, and in which V(t) — e lHot Vs(t)e~ tHot causes transitions between 
eigenstates of Hq- We work in atòmic units, in which e 2 = m e = Ti = 1. 

A. Formulation of the N-body problem in the time domain 

We seek solutions ^f(r%, . . . , r^; t) to the N-body problem described by Eq. Q). It is both conventional and sensible 
to separate the influence of the external interaction V(t) from the initial state at some time to before the external 
fields are applied to the atòmic system, namely, 



*(ri, . ..,r N ;t) = Ü (t, í )*(ri, . ..,r N ; t ) . 



(2) 
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Here all dynamics are contained in the time evolution operator (or Green's function) Ü(t,to). 
It is easily shown that in the intermediate representation U(t, to) satisfies 

ij t Ú(t,t ) = V(t)Ü(t,to). 
The solution for U(t, tg), which may be verified by insertion into Eq. Jjyl, is 



(3) 



U(t,t ) = l-i V(h)dh + {-i) 2 / V(h)dt 



+ 



+ Hf 
(-*)* 



ti 



V(íi)díi / V(t 2 )dt 2 



to 
-i n 



V(t 2 )dt 2 
t n -l 

to 



V(t n ) dt n + 



V{t')dt' 



to 



-Te~ % ^ {t ' )dt ' 



(4) 



Here T is the Dyson time-ordering operator, which arranges the interactions V{t') in order of incrcasing timc, similar 
to the requirements of causality. The time ordering operator is central to the discussion in this paper, as it relates 
both to observable time ordering effects and to time correlation. The key idea is that TV(t\)V(t2) = V(t2)V(tx) if 
Í2 > t\ and V {t\)V (t.2) otherwise. The no time ordering (NTO) approximation is the approximation in which T — ► 1, 
and the constraint of time ordering is not enforced. In this limit all time sequences of the V(t!) are equally weighted, 
so that 



Ü(t,t ) -> C/ N to(í,ío) =e 



-i f* V(t')dt' 

JtQ 



(5) 



in Eq. (gj. 



B. Independent particle approximation (IPA) 

As noted above, the N-body problem is notoriously difficult to solve. A particularly uscful approximation is the 
independent particle approximation (IPA). For an atom with N electrons, H = Ylj(~ V|/2 — Z/rj + J2i<j ^-/ r ij)i 
where Z is the nuclear charge and fj are the electron positions in a coordinate system centered on the nucleus, 
and Vs(t) = ^p/I-^(*) — ^j'I; where Z p and R(t) are the charge and position of the projectile. If ^/ r ij 

is approximated by a mean field y\ v(rj), then H = H a + Vs(t) is reduced to a sum of single-particle terms that 
may be solved using separation of variables. One obtains the independent particle approximation (IPA) to the exact 
N-particle wave function, namely, 

í(r 1 ,...,r J v;í)^n^^' í ) = n^'( í ' í o)^^' í o) ( IPA ) > ( 6 ) 

3 3 

where each ipj(rj,t) is a single-particle wave function. We note here that in the IPA approximation, different times 
may be used for different partides, if so desired, since the partides are independent in both space and time. Whatever 
happens to one of the partides does not influcncc what happens to any of the others, although time ordering is 
retained within the evolution of each independent particle. 



C. Independent time approximation (ITA) 

In order to address the question, "what happens if one tries to define an independent time approximation (ITA), 
where one explores correlations in time rather than the conventional correlations in space?" , a list of comparisons 
between the ITA and IPA was developed 2] . An updated comparison is summarized in Table Q] 

There are similarities between the temporal independent time approximation and the spatial independent particle 
approximation, as seen in Table^ Time and space correlation can each be defined as a deviation from an uncorrelated 
limit, where the uncorrelated limit is given by a product form. Electron identity, which has been ignored here for 
simplicity of presentation, may be restored by antisymmetrizing the uncorrelated single-electron wave functions. The 
uncorrelated limit may also be described by an average of the appropriate correlation operator, as indicated in Table^ 
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TABLE I: Comparison of correlation in space and time. 





Spatial correlation 


Temporal correlation 


Cause: 


spatially varying internal 
Coulomb interactions 


T and V(t) 
time ordering of 
cxtcrnal intcractions 


Origin: 


Ho = Ei Hoj + Y,i<j va 


i£ú = V(t)Ü 


Uncorrelated límit: 


IPA 


ITA 


Product form: 


*(ri,...,r JV ) Uj^( r j) 




No fluctuations: 




T cor = T - T av -> or <JV(í) = í>(í) - V -> 


Average value: 


*" ^av ^mcan íicld 


T — > T av = 1 or V(í) -> V 



Correlation may then be defincd in tcrms of fluctuations away from the average, as is done in statistical mechanics jy] . 
ín both the spatial and temporal cases, the average term may form the basis for useful approximate calculations. 

There are also notable differences Q between temporal and spacial correlation, as detailed in Table While 
correlation in space arises in the asymptotic target Hamiltonian Hq, and affects both the asymptotic initial wave 
function í'(ío) & n d the evolution operator U(t,to), correlation in time occurs only in the time evolution operator 
Ü(t,to)- Correlation in space comes from l/r,j inter-electron interactions within the target. In the IPA, phase 
coherence and time correlation between electrons are both lost, as seen, for example, by noting that matrix elements 

of V in Eq. J7J may be complex and that the time order is significant except when all the V go to a common V . 
Time correlation arises from time ordering of the cxtcrnal intcraction V{t) acting on different partides. The ITA is 
intimately related to the NTO since T — > T av = 1 in both cases. However, in the ITA T — ► 1 is applied only to the 
cross terms affccting different partides, while time ordering for each individual particle is retained. In the ITA, each 
particle evolves independently in time, although the initial state may be spatially correlated. Thus, the initial state 

5Zfc c fcIli ^j^ifjito) evolves to Efe c fcIXjV^ wnere the state of each particle may be evolved using its own 

independent time tj. Removing some or all of the time ordering terms is straightforward in practice since the T cor 
terms are easily identified, at least at second order in perturbation theory 0. 



III. ATÒMIC COL·LISIONS 



In this section we review time ordering and time correlation using perturbation theory for the interaction of charged 
partides with helium 0, • From Eq. Q we see that the leading effect due to time ordering arises at second order 
in the V(t) expansion. Thus the system has an infinite number of states but the external interaction V(t) only acts 
on the system twice. In Sec. IIVI we consider a strongly perturbed qubit, a two-state system that interacts with the 
external field V(t) an infinite number of times. 

Through second order in V(t), the time evolution operator U(t,to) is given by 

Ü(t,t ) = Te" l ^o V(t ' )dt ' ~ l-i í V^dh- í V(ti)dh í 1 V(t 2 )dt 2 

= l-i f V^dh-T^- f V{t x )dh f V(t 2 )dt 2 . (7) 

Jt Z - Jt Jto 

This may also be obtained by integrating Eq. (J3J through second order. Note that as T — > I the order in which 
the interactions act may be interchanged. It is the difference between the limit as T — > I and the full result that 
defines the effects of time ordering. In the case of weak correlation considered here, the NTO limit of T — > 1 will 
coincidentally yield the ITA as well. The trick now is to separate T from T av = 1. 
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A. Primary results 

1. Various pathways to ITA 

To separate the time ordering effects from the non-time ordering (NTO) effects, it is useful to write, 

T = T av + (T - T av ) (8) 

where T av yields the NTO approximation and T — T av yields the effects of time ordering. This decomposition of T 
into an average part plus fluctuations is central for this paper. Note that T av = 1 is required to satisfy the initial 
condition í7(íq,Ío) = 1- 

One conceptual pathway to the ITA proceeds by analogy with the NTO limit. In the ITA, time ordering is enforced 
among all potentials acting on an individual particle (e.g. for Vi{t\)Vi(t 2 ) at second order), but ignored for different 
partides in cross terms such as Vi(t\)Vj(t 2 ) when i ^ j. This distinguishes the ITA from the NTO approximation 
where time ordering is removed for all terms. 

Another clever conceptual pathway to time ordering and the ITA, noticed by Godunov Q , is via use of commutator 
relations. Consider the identity, which we call the Godunov identity [g, 

V(h)V(h) = \ (v(h)V(t2) + V{t2)V{h)) + \[V{h\ V(ía)] ■ (9) 

In the case of atòmic col·lisions with helium, the projectile interacts with both electrons so that V(t') = Vi(t') + V 2 (t'). 
Now [V(ti),V(t 2 )} contains both [Vi(ti),Vi(t 2 )\ and [Vi{h),Vj{t 2 )} {i ^ j) terms. In the NTO approximation, all 
commutator terms are eliminated, i.e. [V{ti), V(t 2 )] -> and V{ti)V(t 2 ) -> \ (v(ti)V(t 2 ) + V'(í 2 )V'(íi)) so that the 
time evolution operator at second order is given by the average value of both possible time sequences. In the ITA, 
only cross commutator terms between fields acting on different partides are neglected, i.e. [Vi(íi), t^- (Í2)] - * and 

V5(ti)V5(í 2 ) -> \ (v t {h)V 3 (t 2 ) + Vji^Mtií) only for i £ j. 

A third pathway to the ITA is through elimination of off-energy-shell effects. The effect of the Dyson time ordering 
operator at second order may be expressed as TV(t\)V(t 2 ) = 0(íi — t 2 )V {t\)V (t 2 ) + (1 «-> 2), and thus 

-t / v{t 1 )v{t 2 )dt 1 dt 2 = / eih-t^vitjvfàdhdh. (10) 

j— 00 j— 00 j— 00 j— 00 

Taking T — > 1 is equivalent to replacing 0(ti — 1 2 ) by the constant 1/2. The Fourier transform of the Heavyside theta 
function is well known, namely 

X) 1 

e «(ti-t a )e( tl _ É2 ) d(íl _ h ) = nS (E) +iP v -. (11) 

> & 

The principal value term accesses off-energy-shell states, which violate energy conservation during the short collision 
time in a manner consistent with the uncertainty relation. Since taking T — > 1 gives the ir5(E) term in Eq. i|ll[l . 
it is the principal value term that carries the effects of time ordering. Ignoring the off-shell contribution (which is 
conveniently n/2 out of phasc with the time-averaged term) whenever the fields V(ti) and V(t 2 ) act on different 
partides yields the ITA. 

There is yet a fourth pathway to the NTO or ITA. For potentials that are time-independent in the intermediate 
representation, [V(íi), ^(£2)] — > and the NTO approximation is exact. In general, the NTO approximation may 
be obtained by replacing the true external potential V(t) with its time average over the duration of the interaction, 
— 1 r* „ 

i.e. V(t') — > V = / V(t') dt' ', so that the exponent J V(t') dt' in the time-ordered exponential of Eq. Q is 

_ t - 1 J to 

replaced by V ■ (t — to). The replacemcnt of the external time-dependent interaction by its time average is analogous 
to the replacemcnt of the true inter-particle interaction by its mean-ficld value in the IPA, which involves averaging 
over the positions of all but one of the partides. 

We note that the limit of constant potential V(t) in which the NTO approximation becomes exact is distinct from 
the adiabatic limit, in which the potential merely changes slowly with time. Furthermore, the content of the NTO 
or ITA approximation depends on the representation used, as we will see explicitly in Sec. IIVI Thus, for a given 
decomposition H = Hq + Vs(t), it is the interaction-representation potential V(t) — é lH ° t Vs{i)e~ lHat that must be 
constant for the NTO or ITA approximation to be exact. In the Heisenberg representation, V(t) — by construction, 
and the NTO or ITA is always trivially exact. To summarize, pathways to the ITA approximation in this section 
include: 
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1. T — > 1, 

2. [K(ti),fy(t a )]-0, 
4. F(í) -> 

The difference between the ITA and NTO approximations is that in the second item above, the commutator disappears 
for all terms in the NTO but only for the cross terms in the ITA. 

2. Economy of NTO 

A variety of second-order calculations with and without time ordering have been done by Godunov and collaborators 
in recent years 0,0] . Most of these calculations have studied two-electron transitions in àtoms caused by high-velocity 
collisions with protons or electrons, where second Born methods are applicable. The advantage of looking at two- 
electron transitions is that second order terms are often dominant, since both elàstic scattering and single-electron 
transitions, corresponding to the first two terms in Eq. J7J, are experimentally eliminated. A disadvantage is that the 
resulting cross sections arc quitc small. The Godunov code is remarkable in that the off-shell terms can be computed 
exactly at second order in pcrturbation theory, in contrast with most othcr cxisting calculations, which use closure 
approximations to avoid this relatively difficult calculation. Dropping the difficult but interesting off-shell terms yields 
a result without any time ordering. Comparing results with and without inclusion of off-shell terms then yields the 
time ordering or time correlation effect. The àlgebra required for the off-shell terms is relatively tedious, as reflected by 
the off-shell calculation typically requiring several hundred times more computer time. Thus the NTO (and possibly 
the ITA) represents a substantial reduction of computer time and algebraic effort. When vàlid, this approximation 
can therefore be used to attack problems harder than those that require the full off-shell terms. Unfortunately, at this 
point no simple, physically transparent criterion is known to us that determine when the NTO or ITA approximation 
is vàlid. 



B. Experimental evidence 

There have been ten or so experiments that show some evidence for time ordering and time correlation effects for 
two-electron transitions in high velocity collisions. Perhaps the most dramàtic evidence is the factor of two difference 
found IlJ in the double ionization of helium in collisions with protons and antiprotons at several MeV. There seems 
to be agreement that this difference is due to time ordering and time correlation, but no definitive theoretical studics 
have yet been done. Some studies of Auger profiles also support the need for time correlation and time ordering. 
The clearest direct comparison of experimental data with time-ordered and non-time-ordered theoretical calculations 
is for polarization of light emitted after excitation-ionization of helium by proton impact |I3| . 

In this study, time ordering or time correlation has a 20% effect on the predicted polarization. The experiment, 
with 5% errors, is in excellent agreement with Godunov's calculations, except at the lowest energy where perturbation 
theory is expected to break down, and clearly shows the importance of time ordering or time correlation effects. There 
is a need for more calculations by different authors to confirm these effects of time correlation and time ordering. 

In all of these two-electron transition experiments, one may do a simultaneous expansion in the external interac- 
tion V(t) and the correlation interaction v. Time correlation effects are associated with commutators of the form 
[Vi(íi), Vjfe)], where i and j label two different electrons, and are second order in V. In addition, there are time 
ordering effects not related to time correlation: these are associated with commutators [V^(íi), Vi(Í2)] and enter at 
second order in V and first order in v (since the commutator vanishes for v = 0). Thus the difference between time 
ordering and time correlation is a third order effect, smaller by at least a factor of 10 than either the time ordering 
effect or the time correlation effect taken individually. The experiments show direct evidence for time correlation but 
do not distinguish between time correlation and time ordering. 

There is also one clear and definitive study by the group of Thomas 0] that shows direct evidence for time ordering 
in àtoms interacting with a time-varying magnètic field. 



8 



0.25 
0.20 



0.15 

| 0.10 

to 

£ 0.05 
_ço 

°- 0.00 
-0.05 
0.10 

2 3 4 5 6 7 8 
collision velocity (a.u) 

FIG. 1: Calculations with and without time ordering between electrons are compared to experimental data. Here polarized light is 
emitted from helium following ls — > 2p excitation of one electron accompanied by ionization of the second electron ll.'.fl . The polarization 
fraction is plotted as a function of the velocity of the incident proton. The first-order calculation (Born 1) has no time ordering. The 
second-order calculation is shown both with time ordering (Born 2 full) and in the no time ordering (NTO) approximation (Born 2 unc). As 
explained in the text, the NTO approximation is the same as the independent time approximation (ITA) at this order since the interaction 
is weak. 



IV. QUBITS 

A qubit is a very simple two-state (e.g. on and off) quantum system whose state population may be changed 
by an external potential V(t), in analogy with the way in which an atòmic state may be changed by the Coulomb 
potential V(t) of the projectile in an atòmic collision. Qubits are building blocks for the complex interconnected 
N-qubit systems that can be used for quantum computation and quantum information. The extension from a singlo 
qubit to a system of N interconnected qubits is analogous to the extension from a one-electron atom to a correlated 
N-electron atom. The interaction of qubits with each other and with their environment still have to be dealt with 
before a quantum computer becomes a reality. 

In this section we consider time ordering in a single qubit, analogous to time ordering in scattering from atòmic 
hydrogen. To better understand time ordering effects, we work in the time domain rather than the more commou 
energy (or frequency) domain. The advantages of working with qubits include the possibility of easily handling non- 
perturbative external potentials and the existence in some cases of analytic solutions, so that numerical calculations 
can be avoided. This yields new ways to think about time ordering mathematically and physically. Specific effects 
due to time ordering in a simply pulsed qubit are shown, for example, in Fig.[3]below. The idea is to extend the above 
analysis of time ordering for weakly perturbed atòmic collisions to the case of strongly perturbed qubits. In order 
to make use of transparent analytic expressions, we spend some effort discussing the time evolution of qubits pulsed 
sharply in time, commonly referred to as "kicked" qubits. With the exception of the last brief subsection, there is 
little emphasis here on time correlation between qubits, simply because little work has been done on this problem 
to our knowledge. So our emphasis in this section is almost entirely on time ordering. One new feature here will 
be to explore what happens to time ordering under change of representation - specifically when we change from the 
Schròdinger to the intcrmediate representation. 

A. Overview of single qubits 

The qubit wave function ijj(t) is a linear superposition of the "on" and "off" states, namely 
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a 2 {t) 



(12) 



Population can be transferred from one state to the other by applying an external potential V(t), which can have 
the form of a single pulse characterized by a time duration, r. The full Hamiltonian can be written in terms of the 
Pauli spin matrices, 



H{t) = H + V s (t) 



AE/2 


+ 


V{t) 


+AE/2 




V(t) 



= — (Jz + V{t)a x . 

The time dependence of the system can be described by the time evolution operator U(t,to). 
simplify the notation, we can write tp(t) = U(t)ip(0), where 



(13) 



Allowing to = to 



U(t) = 



U n (t) U 12 (t) 

Ü2l{t) U 22 (t) 



= Te 



-i J* V(t')dt' 



(14) 



Here V{t')= e lHat V s (f 



,-iiïní' 



is the interaction potential in the intermediate representation, and U(t) is the formal 



solution to the differential equation In the Schròdinger representation, V(t') is replaced by H(t') = Hq + Vg(í')- 
The time evolution of a qubit depends on the energy splitting AE and the time dependence V(t) of the external 
potential. Depending on the complexity of V(t), U{t) may or may not have a simple analytic form. In the latter case, 
numerical solutions of the coupled differential equations 10 may obscure information about these quantum systems. 
Analytic solutions, when available, are more convenient and easy to analyze. 



B. Simply pulsed qubits 

Here we consider simply pulsed qubits - that is, qubits subject to an external field V{t) that has a finite duration 
in time and a sensibly simple shape, such as a simple rectangular or Gaussian pulse. Such pulses are convenient for 
studying qubits in the time domain. In some cases they also lead to convenient analytic solutions, even for strong 
fields. 



1. Qubit map 



A qubit map, such as that shown in Fig. [3 is a tool enabling one to visualize how the possible behavior of simply 
pulsed qubits depends on the variables AE and V(t). In order to make the qubit map a helpful visual tool, the 
map coordinates are taken to be a dimensionless level splitting Ai? r/2 (where r is the time duration of the pulse 
V{t)) and a dimensionless pulse strength V(t') dt' . These two coordinates determine the effect of the unperturbed 
Hamiltonian Hq and of the external potential V(t) on the evolution operator U(t). It is useful to think of these two 
variables as independent phase angles or action-like integrals. 



2. RWA solutions 



The kicked qubits we focus on below provide an alternative to the well-established rotating wave approximation 
(RWA) method |15l lla Il7| based on a single resonant transition frequency, that may be detuned. Our emphasis 
here is on kicked (i.e. strongly time-localized) potentials because: i) the physics is then naturally analyzed in the 
time domain, and ii) kicks have been less widely explored than the RWA approach, which works well for sharp pulses 
in the reciprocal frequency space. While the RWA method is useful for two-state systems perturbed by an external 
interaction of narrow bandwidth, it fails to describe simply kicked two-state systems where the pulse bandwidth may 
be very broad. 
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FIG. 2: Qubit map for qubits interacting with simply pulsed finite external potentials. Here r is the duration of the pulse, AE is the 
energy difFerence between the two states of the qubit, and V(t) is the external potential. When the total phase associated with the external 
potential is small, i.e. V(t')dt' -C 2?r, then the expression for the time evolution operator may be expanded in powers of V using 

e —iV(t )dt i _ iV(t')dt' and only the first few terms retained. This corresponds to Standard perturbation theory for either quickly 
perturbed qubits where r is small or slowly perturbed qubits where r is large. Similarly, if the phase associated with the energy splitting 
is small, AEt/2 -C 2tt, then we can treat the qubit as degenerate, and the solution can be expanded in powers of AEt. If both phases 
are large, then the adiabatic approximation generally applies. 



3. Kicked qubits 



A useful approximation for simply pulsed qubits is the fast, narrow pulse or "kick" limit in which the width r of 
the pulse goes to zero, while the integrated strength or area under the external potential curve 

t 

V(t')dt' (15) 

remains fixed. Formally, the shape of a very narrow pulse of finite total strength a may be expressed by a delta 
function: V{t') = aS(t' — tk), where tk is the time at which the pulse is centered. The kicked region corresponds to 
the lower half of the qubit map in Fig.|2] Here the duration of the pulse is so short that Ai? r/2 <C 2tt, i.e. there is 
not enough time for the splitting AE to have a significant effect while the pulse is active. The integrated strength of 
the pulse, a, may be either large or small in this region. If a is large, we are in the lower right quadrant of the map, 
where the kicked region overlaps with the adiabatic region. If a is small, we are in the lower left quadrant, where the 
kicked region overlaps with the perturbative region. 

Single kick 




For a two-state system subjected to a single kick at time tk, corresponding to V(t') = aS(t' — tk), the integration 
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over time is trivial and the time evolution operator in Eq. I|14[) becomes 



Ü k (t) = Texp 

= exp 

= exp 

= exp 



■ / e - l AE^t'/2 a§ ^ _ tk)(Jxel AEa z t'/2 dtl 



r iAEt t /2 









1 



-*Q 1 e iAEt k 



cos a 

iAEt 



gíAiï t k /2 ) [l Q 
iAEt k 


-iAEt, 



iAEt k /2 







-iAEt k /2 



sin a cos a 



(16) 



for í > ífc. The last line of Eq. (jl6|l was obtained by expanding the fourth line in powers of a using the identity 



-iAEt k 



e 

iAEt k q 



2ii 



= /. Equivalently, one can take advantage of the useful identity (which we also used in the 



third line) e 1 ^ <T " = cos cf> ï + i sin <fi a ■ u, where ü is an arbitrary unit vector. Note that Ü k (t) is independent of the 
final time t in the intermediate representation. 

The occupation probabilities for a kicked qubit initially in state 1 are 



Pi (i) = | ai (í)| 2 = l^íi(í)| 2 = cos 2 a 

P 2 (t) = \a 2 (t)\ 2 = |t4(t)| 2 = sin 2 a . (17) 

This simple example may be extended to a series of kicks 18] . It is one of the few cases in which analytic solutions 
may be obtained for qubits controlled by external potentials. 

Double kicks 

The simplest example of a series of arbitrary kicks is a sequence of two kicks of strengths ai and a 2 , applied at 
times ti and t 2 respectively, i.e. Vs(t) = (a\8(t — t\) + a 2 S(t — t 2 ))a x . Eq. l|14fl is then easily solved in the interaction 
representation, namely, 

fjkzM = jjk 2 x fjk, 

_ í cos a 2 
~ \~ie lAEt2 sma 2 

( U 1X U 12 \ 

\U 21 U 22 ) ' 

where 

Un = cosai cosa2 — sin ai sina2 e ~ lAEt - ; (19) 

U 21 = - le iA£t +/ 2 (cosa 1 sina 2 e ïAí;t -/ 2 +sina 1 cosa 2 e- ïAí;t -/ 2 ) . 

Here í_ = t 2 — íi, and t+ = t\ + t 2 . In the limit t 2 — »• t%, Eq. I|18|l reduces to Eq. Hlfi|l with a — > ol\ + a 2 . Note that 
[U k2 , U kl ] ^ so that the time ordering of the interactions is important. 

The àlgebra for a combination of two arbitrary kicks |19| , one proportional to a y and the other proportional to a x , 
is very similar to the above. Triple kicks are also straightforward to solve analytically. 



_ ie -iAEt 2 ginQ , 2 

cos a 2 



cos «i 
-ie lAEtl sin ai 



_ ie -iAEt! s [ nai 

cos ai 



(18) 



4- Time ordering in a doubly kicked qubit 



In this subsection we use our analytic expressions to examine the effect of the Dyson time ordering operator T in a 
kicked two-state system. Time ordering has been considered previously in the context of atòmic collisions with charged 
partides 0, U 0[ ÍÏH |20| and differs somewhat from the order in which external puls es are applied, as illustrated 
below. As is intuitively evident, there is no time ordering in a singly kicked qubit |2l| since there is only one kick. 
The simplest kicked two-state system that shows an effect due to time ordering is the qubit kicked by two equal and 
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opposite pulses labeled k and — k separated by a time interval í_ = Í2 — íi. The evolution matrix for this system 

is Un, 

~ /e- ^A£t -/ 2 (cos4^_+^cos2asin^i_) e - <ABt + sin2asin 4^í_ \ 

( ~ ^ -e lA£t +sin2asin^í_ e lAE '-/ 2 (cos - i cos2a sin ^í_) J ' ( j 



The time evolution operator ÍA°) in the limit of no time ordering, i.e. in the approximation T — ► 1, may in principio 

be generally obtained by replacing f* V(t') dt' with Vt, where V is an average (constant) value of the interaction. 
In our case, it is then straightforward to show 



T y<o)-k,k _ -iVt _ cos^asm— 1_) e — "+ sm{Za sm — í_ j \ 

U ~ e ~\ „iAEt+ e - fr, AE . \ „„„Co„, A_B . \ | ■ V zí ) 



cos(2a sin 4^í_) e~ lAiï '+ sin(2a sin ^-tJ) 

_ e iAEt+ S i n (2asin ^í_) cos(2asin ^í_) 

In this example we now have analytic expressions for the matrix elements of both the evolution operator (j~ k ' k that 
contains time ordering and the evolution operator JJ í S ) )^ k ^ k without time ordering. 




FIG. 3: Difference in population transfer probability, P2 — P 2 vs - e = sin(A_E Í-/2) and <fi = 2a. Here t_ = Í2 — íi is the time between 
the pulses and a = J V(t')dt' is a measure of the interaction strength. The two-state system is kicked by a sharp pulse of strength a at 
time íi and by an equal and opposite pulse at time Í2- The difference, P2 — is due to time ordering in this qubit. 



It may be shown analytically |ljj that for two kicks both proportional to a x , the order of the kicks does not change 
the final transfer probability Pi. However, interestingly, this does not mean that there is no effect due to time ordering 
in this case. As we show next, there is an effect due to time ordering in this case, even though interchanging the order 
of the kicks has no effect. The effect of time ordering on the occupation probabilities may be examined by considering 
the probability of population transfer from the on state to the off state with and without time ordering, namely from 
Eqs. (23 and 

AE 

P2 = \U 2 i\ 2 = | sin 2a sin— í_| 2 = |esin0| 2 , (22) 
= |^? ) | a = |sm(2a 8 in^i_)| 2 = |sinc^ a , 

where e = sin 4pí- and <f> — 2a. 

The effect of time ordering is shown in Fig.|3J where P2 — P 2 1S plotted as a function of (f> = 2a, corresponding to 
the strength of the kicks, and e = sin =¥-t-, which varies with the time separation of the two kicks. The effect of time 
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ordering disappears in our example in the limit that either the interaction strength or the time separation between 
the pulses goes to zero. For small, but finite, vàlues of both the interaction strength and the time separation between 
the pulses, the effect of time ordering is to reduce the transition probability from the initially occupied state to an 
initially unoccupied state. That is, in this regime time ordering reduces the maximum transfer of population from 
one state to another. As either of these two parameters gets sufficiently large, the effect of time ordering oscillates 
with increasing vàlues of the interaction strength or the inter-pulse separation time. Time ordering effects are present 
even though U~ k > k = JJ k ~ k . 

5. Time ordering vs. time reversal 

Let us now pause to examine the difference between time ordering and time reversal in this simple, illustrative 
example. Reversal of time ordering means that, since kick strengths otk and kick times tk are both intcrchanged, 
í_ — * — t- and a — > —a. In this case one sees from Eqs. H20|) and Ij21ll that [/~ fe > fe experiences a change of phase, 
while Ü^~ k ' k is invariant under change of time ordering. For time reversal [22^. t± — ► —t± and, since the initial and 
final states are also interchanged, U — > W . Inspection of the same equations as above shows that TJ~ k ' k and U^~ k ' k 
are both invariant under time reversal, as expected. When the symmetry between the kicks k\ and ki is broken, i.c. 
a 2 ^ ±a x , the difference between Ü k2 ' kl = Ü k2 Ü kl and Ü kl ' k2 = Ü kl Ü k2 can be observed [ïjj. 

6. Numèrica! calculations of time ordering 

As an illustrative specific example, we present the results of numerical calculations for 2s — > 2p transitions in atòmic 
hydrogen caused by a Gaussian pulse of finite width r. The occupation probabilities of the 2s and 2p states are 
evaluated by integrating the two-state equations using a Standard fourth order Runge-Kutta method. This enables us 
to verify the validity of our analytic solutions for kicked qubits in the limit r — > and also to consider the effects of finite 
pulse width. In this system, the unperturbed level splitting is the 2s — 2p shift, Ai? = E 2p ~E 2s = 4.37 x 10~ 6 eV. The 
corresponding time scale is the Rabi time, Tae = 2ir/AE — 972 ps, which gives the period of oscillation between the 
states. In our numerical calculations we use for convenience a Gaussian pulse of the form V(t) = (a/y / 7rr)e _< · t_í ' e - ) l T . 
The two-state coupled equations for the amplitudes a± and a 2 of Eq. 112|) take the form implied by Eq. (|13|) . 

ià x = -lA£ fll + -^ (í ~ tfc)2/r2 a 2 

2 ^HT 

ià 2 = \AEa 2 + ^e- {t - tk)2/T2 ai . (23) 

Here the pulse is applied at tf. = 150 ps and we have chosen a — n/2 so that in the limit of a perfect kick all of the 
population will be transferred from the 2s to the 2p state after t = í/.. In Fig.^J the resulting transition probability 
is shown as a function of pulse width r and as a function of observation time Tf for Tf > í& . 

In the Schròdinger picture there are very large differences between the transition probabilities with and without 
time ordering, P 2 (T f ) and p!f\Tf), even for an ideal kick. This occurs because the energy splitting AE is non-zero, 
and for Tf > a/AE — cíT&eI'Z'k, the average potential V = a/Tf becomes smaller than the energy splitting AE. 
Thus, for a given pulse, the influence of the potential necessarily decreases at large Tf, and any transfer probability 
becomes exponentially small. In effect, the free propagation before and after the pulse diminishes the effect of the 
pulse itself in the Schròdinger picture, when time ordering is removed. This behavior contrasts with the intermediate 
picture result, where Pf 2 {Tf) depends on AEt/2 but not on Tf, as seen on the left side of Fig. 0] The contrast is 
also evident on the right hand side of Fig. 21 where, after the pulse has died off, the value of P 2 °\Tf) decays with 
increasing Tf, while Pj° 2 \Tf) approaches a constant. 

C. Networks of qubits 

To our knowledge relatively little has been understood analytically for systems of coupled qubits. Understanding 
coupling between qubits is a well recognized challenge in the field of quantum computing |23j. In our view, developing 
a realistic analytic model for two coupled qubits could provide a useful and instructive example in the fields of 
quantum computing, quantum information and coherent control. Unitary evolution operators acting on a system of 
non-interacting qubits formally belong to the SU(2) x SU(2) x • • ■ x SU(2) unitary group, which is simply the set 



14 




FIG. 4: Target state probability as a function of the pulse width r (on the left) and of observation time Tf (on the right). Here 
Tab = 2tt/AE is the Rabi time for oscillations between the states, where AE = E2 P — E^a- The heavy line denotes probability including 
time ordering, the dashed line denotes the probability in the Schròdinger picture without time ordering, and the dotted line represents the 
probability in the intermediate picture without time ordering. On the right, the lines begin at the midpoint of the pulse, Te = tfe. The 
Schròdinger results damp out for large Tf as explained in the text. 



of all local qubit operations. This set is a subgroup of the full SU(2 N ) dynamic group of N coupled qubits. The 
SU (4) dynamic group of two interacting qubits plays a fundamental role in the analysis of multi-qubit dynamics since 
any operator (i.e. quantum gate) from the full SU(2 N ) group can be factorized as a product of SU(4) two-qubit 
gates. In connection with this property, it is useful to note that any 4-level quantum system can be used to encode 
a tensor-product four-dimensional Hilbert space of two qubits. The specific form of such encoding is completely 
determined [24| by fixing one point on the orbit of the maximal SU(2) x SU(2) subgroup of SU(A). Next, the 3-qubit 
system (e.g. a carrier space for GHZ multiparticle entangled states) can be similarly encoded in an arbitrary 8-level 
quantum system. This requires two steps: identifying the maximal SU(4) x SU{2) subgroup of SU(8), and then 
adjusting the SU(2) x SU(2) subgroup of the resulting SU{4) group. 

It would be useful to have a well developed theory for time correlations between interacting qubits - one that 
clarifies how the time dependence of a field acting on one qubit impacts the time evolution of another qubit, for 
example in switching. However, defining the independent time approximation for N-qubit systems in a useful way is 
a challenge. Although working with n = 2 N degenerate states is doable in principle [2!| (setting aside the problem of 
solving an n-th order equation for n > 4) , this does not always give the NTO approximation, seemingly a prerequisite 
for the ITA approximation needed to define time correlations. 

Furthcrmore, in many cases sequencing of external interactions can be problemàtic. For example, if V(t') = 

Va(í') + V B (f) for two partides or qubits A and B, and Ü = Te^ l I {VA{t ' )+VB{t ' ))dt ' , then the ITA is not Ü A ■ Ü B in 
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general. The ITA is given by Üa ■ Üb = Üb ■ Üa if [Va(í), Vb(í')] = 0, but in that case all time correlations vanish and 
the ITA is exact. Also, it difficult for us to envision how one may satisfy [Ua{í), ?/b(í')] = with [Uj(t), Uj(t')] ^ 
for J = A, B. That is, in what situations can one eliminate inter-particle time correlations while retaining time 
ordering for individual partides? If all commutators terms vanish, then time correlation cffects and time ordering 
for individual partides are both absent. Similarly, time correlation effects disappear if all parts of a composite V(t) 

are replaced by the time averaged value V, but this again is equivalent to eliminating all time ordering, even within 
single-particle evolution. 

V. SUMMARY 

In solutions of the time dependent Schròdinger equation there are only two sources of time dependence, namely the 
Dyson time ordering operator T and the explicit time dependence of the interaction V(t). This simplifies the study 
of how time dependence may influence the evolution of A^-body quantum systems. The causal-like constraint of time 
ordering between fields acting on different partides can cause time correlation between the partides. That is, the time 
dependence of a field acting on one particlc can influence the evolution of other partides. Correlation is traditionally 
studied by defining an uncorrelated limit. In the case of time correlation, we have called this the independent time 
approximation (ITA), and have pointed out similarities (e.g. in Table[ÏJ to the widely used and practical independent 
particle approximation (IPA) that eliminates spatial correlations between partides. Similarly the limit of no time 
ordering (NTO) can be defined by eliminating all time ordering constraints, for fields acting on the same particle or 
on different partides. Thus, the ITA may be viewed as the NTO applied to cross terms only. 

The ITA or NTO limit may be reached in several ways, but the most general seems to be to define a mean time- 
averaged coupling interaction (as is done in the IPA). Time ordering effects have been observed in weakly perturbed 
atòmic collisions. One may also consider strongly perturbed systems of coupled qubits. However, relatively little has 
been done on this problem, which is a key problem in quantum computing. In this paper we considered the effect of 
time ordering on strongly perturbed single qubits. To do this, we focused on qubits subjected to fast strong external 
pulses called kicks, where useful analytic expressions for observable transition probabilities may be obtained. What we 
have found (but not discussed here) is that the NTO may be easier to implement than the ITA. This means that the 
ITA might be computationally awkward in general. We have also demonstrated that the NTO (and consequently the 
ITA) are dependent on the representation used. There is evidence that the intermediate representation is preferred. 
This suggests that one could find gauge dependence in specific NTO and ITA terms if MBPT is used, and raises a 
question 0, |2(j about the physical meaning of time ordering and time correlation. 

In summary, the methods we have developed in this paper are intended to probe the nature of how time works 
in quantum N-body systems. The focal point in our approach is the influence of the constraint imposed by time 
ordering, which can lead to time correlation between different parts of the system. 
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